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Quantum metrology for the Ising Hamiltonian with transverse magnetic field 

Michael Skotiniotis, Pavel Sekatski, and Wolfgang Diir 
Institut fiir Theoretische Physik, Universitdt Innsbruck, Technikerstr. 25, A-6020 Innsbruck, Austria 

(Dated: August 26, 2015) 

We consider quantum metrology for unitary evolutions generated by parameter-dependent Hamil¬ 
tonians. We focus on the unitary evolutions generated by the Ising Hamiltonian that describes the 
dynamics of a one-dimensional chain of spins with nearest-neighbour interactions and in the pres¬ 
ence of a global, transverse, magnetic field. We analytically solve the problem and show that the 
precision with which one can estimate the magnetic field (interaction strength) given one knows the 
interaction strength (magnetic field) scales at the Heisenberg limit, and can be achieved by a linear 
superposition of the vacuum and N free fermion states. In addition, we show that GHZ-type states 
exhibit Heisenberg scaling in precision throughout the entire regime of parameters. Moreover, we 
numerically observe that the optimal precision using a product input state scales at the standard 
quantum limit. 

PACS numbers: 03.67.-a, 06.20-f, 75.10.Pq 


I. INTRODUCTION 

Quantum metrology is one of the archetypical applica¬ 
tions where quantum mechanics demonstrably exhibits a 
vast improvement over the best known classical strate¬ 
gies. The use of entangled input states of N qubits, such 
as GHZ states [T], is known to allow for a precise deter¬ 
mination of an unknown parameter such as the relative 
phase in a Mach-Zender interferometer mg, or the fre¬ 
quency of an atomic transition mm with a precision 
that scales inversely proportional to N, the Heisenberg 
Jimit. In comparison the best classical estimation strat¬ 
egy, which employs separable states, gives a precision 
scaling inversely proportional to \/lV, the standard quan¬ 
tum limit [ZlIH]. This observation has proven very useful 
in the development of ultra-precise atomic clocks HE], 
high resolution imaging [T2|, as well as the detection of 
gravitational waves [131 Ej • 

In the absence of any noise or decoherence effects the 
parameter of interest. A, in a metrological scenario is im¬ 
printed onto the state of N probes via the unitary op¬ 
erator U = where H{\) is the Hamiltonian de¬ 

scribing the dynamics of the N probes. For the case of 
phase and frequency estimation—where the parameter to 
be estimated is Xt and A respectively—the parameter is 
a multicplicative factor of the Hamiltonian, H{X) = XH, 
and the Hamiltonian is local, i.e., H = where 

is the Hamiltonian describing the evolution of the I*** 
probe. Indeed, quantum metrology using local Hamilto¬ 
nians where the parameter of interest enters only as a 
multiplicative factor have been studied extensively both 
in the absence as well as in the presence of noise [T5IE5] . 

In stark contrast quantum metrology with more gen¬ 
eral Hamiltonians is only now beginning to attract atten¬ 
tion. Some instances of quantum metrology with param¬ 
eter dependent Hamiltonians concern the estimation of 
time-varying signals [Ml - ES] , the estimation of magnetic- 
field gradients along a spin chain [27ll29j . or the estima¬ 
tion of the anisotropy and/or decoherence along a spin- 
chain using nonequilibrium states |30] . The general prob¬ 


lem of performing quantum metrology using parameter 
dependent Hamiltonians was treated in iniisij where it 
was shown that for parameter dependent local Hamil¬ 
tonians Heisenberg scaling precision with respect to the 
number of probing systems is possible. 

In this work we focus on noiseless quantum metrology 
of N interacting probe systems described by the Ising 
Hamiltonian 

N-l N 

H{J, B) = jY, ^ (I) 

i=l i=l 

where we are interested in determining the precision with 
which one can estimate either the strength of the trans¬ 
verse magnetic field, B, or the coupling interaction, J, 
provided the remaining quantity is known. The Hamil¬ 
tonian of Eq. 0 is known to exhibit a phase transi¬ 
tion |33|, which have been discussed previously in the 
context of quantum metrology and were shown to be re¬ 
sourceful |34l-[37|. Moreover, as the Ising Hamiltonian is 
entanglement generating |3H|, it has found applications 
in ion-trap quantum computing architectures [sa io], 
where either J or B can be controlled at will by modifying 
either the separation of the ions, or the global magnetic 
field strength. 

Here we report the following results regarding precise 
estimation of parameters of the Ising Hamiltonian: 

1. The ultimate precision with which one can estimate 
either J or B, having complete knowledge of B or 
J, using N probe systems scales at the Heisenberg 
limit, and we provide an analytic expression for this 
achievable precision as well as the optimal states 
that achieve it. Our result answers the conjecture 
in |3Ij that Heisenberg scaling is the ultimate limit 
also for the case of N interacting probes for the 
case of nearest-neighbour interactions. 

2. We provide numerical evidence that the ultimate 
achievable precision strictly outperforms the opti¬ 
mal classical strategy, which deploys the N initial 


2 


probes in a pure product state. For up to = 11 bound m 
our numerical study shows that the optimal prod¬ 
uct input state yields a precision that scales at the 
standard quantum limit. 


6Xi > 


1 


( 2 ) 


3. We analytically derive the precision achieved by the 
GHZ-type states, known to achieve the optimal pre¬ 
cision when either J or B are equal to zero, and 
show that these states retain their Heisenberg scal¬ 
ing in precision (up to a constant factor) over the 
entire regime of parameters. 


where J^{px,t) is the quantum Fisher information (QFI) 
of the state pA,t describing the N probes after the unitary 
dynamics have acted. In the most general case the QFI 
can be computed as [41] 

•^(PA,t) = Tr{LAiPA,tAAj , (3) 


This paper is organised as follows. In Sec. [IT] we re¬ 
view the basics of quantum metrology as well as some 
important mathematical results regarding unitary oper¬ 
ators generated by parameter dependent Hamiltonians. 
In Sec. |HI| we use the Jordan-Wigner transformation to 
determine the maximal possible precision with which one 
can estimate either J and B using N systems as probes. 
We compare this to the best possible precision that can 
be achieved by a separable state of N probes which we de¬ 


termine numerically for up to = 11 qubits. In Sec. IV 


we analytically determine the performance of states that 
are known to be optimal at the extreme cases where 
J = 0 and B = 0. Finally, Sec. |V| contains the conclu¬ 
sions of our investigation as well as some open questions 
for future work. 


II. BASICS OF QUANTUM METROLOGY 

In this section we review the main results in quantum 
metrology and review some important facts pertaining to 
parameter-dependent Hamiltonians in general and to the 
Ising Hamiltonian (Eq. Q) in particular. Specifically, 
we concentrate on noiseless quantum metrology and the 
quantum Fisher information (QFI); the central quantity 
of interest in quantum metrology. After introducing the 
QFI we provide a formula for calculating it for the case of 
general parameter-dependent Hamiltonians and then to 
the specific case of the Ising Hamiltonian. We note that 
the behaviour of the QFI, and in particular its scaling 
with time and number of probe systems, was investigated 
for general parameter dependent Hamiltonians in [311132] . 

A standard protocol in noiseless quantum metrology 
can be formulated as follows: N probes are prepared in 
a suitable state and undergo an evolution for some time, 
t, described by the unitary operator t7(A,t) = 
where H (A) is the Hamiltonian describing the dynam¬ 
ics of the N probes and explicitly depends on the vec¬ 
tor of parameters A = (Ai,...,Am)- Finally the N 
probes are measured and an estimate of A is obtained 
from the measurement statistics of v repetitions of the 
above procedure. In what follows we shall assume that 
all other parameters except A^ are known, and shall be 
concerned with estimating Ai as precisely as possible. A 
lower bound on the error, 6Xi = Xi — Xi, for any unbi¬ 
ased estimator Xi is given by the quantum Cramer-Rao 


with 


q:,/3 

the symmetric logarithmic derivative, where a (/3) are the 
eigenvalues of px,t, IV'a), dV'/?)) the corresponding eigen¬ 
vectors, and the sum in Eq. Q is over all a, f3 satisfying 
a + P ^ 0. Here and in what follows dx = -§^- We now 
review the case where the parameter of interest enters as 
a multiplicative factor of the Hamiltonian 


A. Parameter independent Hamiltonians 


In the case of noiseless metrology, an easier expres¬ 
sion for computing the QEI exists if one initializes the N 
probes in a pure state p = In this case the QEI 

can be shown to be 

fF (|^A,t)) = 4 [{^Xi'^l)X,t\^>‘^i’>^■.^) - l(9Ai'0A.t|'0A,t)|^) , 

(5) 

where |'0A,t) = U{\,t)\ijj). If H{\) = XiH then a bit 
of algebra yields fF {\'4’x,t)) = where = 

\ip) — H \il))\^ is the variance of H with re¬ 
spect to the state |i/;). If, in addition, H = X]i=i 
where Vi, j is the Hamiltonian act¬ 

ing on the probe system, then it can be shown that 
when l-ip) = are the 

eigenstates corresponding to the minimum and maximum 
eigenvalues of h, then F (IV’A.t)) = i^-^(amax - amin)^ 
and give the standard quantum limit in estimation pre¬ 
cision |5|. On the other hand if the probes are prepared 
in the Greenberger-Horne-Zeilinger (GHZ) state 



(|an 


1 Q!n 


. . . , Q^max) ) ; (6) 


N times 


N times 


then F{\tpx,t)) = - a^inY, the Heisenberg 

scaling in estimation precision [Sj. In general, Heisen¬ 
berg scaling in precision implies an improvement in scal¬ 
ing with regards to the number of probe systems, over 
the optimally achievable precision which uses the same 
number of probe systems initialized in a separable state. 
Hence, in the case where the parameter to be estimated 
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is a multiplicative factor of a local Hamiltonian, the use 
of highly entangled states leads to a quadratic improve¬ 
ment in scaling precision. Note that besides the GHZ 
states there exists a large class of pretty good states that 
scale at the Heisenberg limit up to a multiplicative fac¬ 
tor [23]. 


B. Parameter dependent Hamiltonians 


We now use Eq. ([^ to compute the QFI for Hamilto¬ 
nians of the form 


H{X) = AiETi + X2H2 ( 7 ) 

where [Hi,H 2 ] 7 ^ 0. We make no assumptions on the 
structure of Hi, H 2 ] in particular we do not assume 
that they are local Hamiltonians. Notice that the Ising 
Hamiltonian of Eq. ([T) is a special case of Eq. 0 . 
whose QFI was studied in detail in m with XiHi — 

J and X 2 H 2 = B YliLi ■ From Eq. 0 

we need to compute \d\^ip\^t) = d\.U{X,t) Itp). As 
[Hi,H 2 ] 7 ^ 0 and using 


d(l-^^iXiHi + X2H2))"' 

- - = lim —^- - ---^— 

dx N^oo Ox- 

(* 

Eq. reads 


^ (IV'A.t)) = 4 ((V^l C/t (A, t)Of{X, t)UiX, t) \^) 
-\{^\U\X,t)0,{X,t)U{X,t)\i;)\") 

= AA^O,{X,t), (9) 


where 


0,{X,t) 



ds U{X, s)HiU\X, s) 


( 10 ) 


III. ESTIMATION OF MAGNETIC FIELD AND 
INTERACTION STRENGTH 

We are interested in determining the optimal precision 
in estimating either the magnetic field strength, B, or in¬ 
teraction strength, J, of the Ising Hamiltonian (Eq. 0 ). 
In particular, we will show that the optimal precision in 
estimating either B or J scales at the Heisenberg limit, 
up to a constant factor which depends only on the ratio of 
J and B, and is achievable by states that are linear super¬ 
positions of the vacuum and fully occupied states of free 
fermions of a suitable type. Furthermore, we numerically 
determine the best achievable precision using a separa¬ 
ble state for up to 11 qubits and show that the optimal 
precision scales, to within best ht errors, linearly with 
N. Hence, our results provide strong evidence that the 
entanglement generated by the Ising Hamiltonian when 
acting on an initially pure separable state is not enough 
to boost the precision in estimation from the SQL to the 
Heisenberg limit. 

We begin by first determining the optimal precision 
in estimation of either B, or J, and the corresponding 
optimal states. To do so we note that via the use of the 
Jordan-Wigner transformation |?2] the Ising Hamiltonian 
of Eq. 0 can be expressed as a quadratic Hamiltonian in 
fermionic creation and annihilation operators which can 
be suitably diagonalized. Specifically the mapping 



and its inverse 



and the variance of Oi is computed with respect to the 
evolved state IV’x.t)- 

Thus, the QFI is maximised by the states that are 
linear superpositions of the eigenstates corresponding 
to the minimum and maximum eigenvalues of Oi{X,t). 
An interesting question is whether the presence of AiiLi 
{X 2 H 2 ) can help boost the precision of estimation of pa¬ 
rameters for A 2 (Ai) respectively. It was shown in [31] 
that this is not the case. Specifically, if one has control 
over either AiiJi or X 2 H 2 and wishes to estimate A 2 or 
Ai respectively, then the optimal strategy is to set the 
dynamics under our control to zero. 

In the next section we use the expressions in Eq. 0 to 
determine the optimal QFI for either the magnetic field B 
or interaction strength J of the Ising Hamiltonian using 
separable and entangled states respectively. 


where {cx, <Jy, tr^} are the Pauli matrices [43] . with 
= ax'^ ± icry\ and aj, are the fermionic anni¬ 
hilation and creation operators for mode j respectively 
and satisfy the anti-commutation relations {a^ oj,} = 


{ajQk} = 0 , {aj, 4 } 

Eq. 0 yields 


Sjkt- Substituting Eq. (12) into 


N N 

H{J,B) = J ^(4 - + aj+i) + 2B'^a]aj. 

i=i i=i 

(13) 

Any quadratic Hamiltonian in the fermionic operators 
can be brought to the diagonal form 


N-l 


H{J,B) = 2^yjal + Pi clck, 


(14) 
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where ttk = Jcos(^) + B, f3k = Jsin(^), and 4, 
Cfc are the creation and annihilation operators of free 
fermions in mode k, with c^r = Cq. The eigenstates of 
H{J,B) are fermionic Fock states, \k), where k is an N- 
bit binary string indicating which modes are occupied 
by fermions. Without loss of genera lity we m ay set the 
vacuum state, |0) to have energy a/«o + <®o = 0- The 
maximally occupied fermionic Fock state, |1) has energy 

equal to Ylh-^ n \/^l+ Pi/ 

As Eq. ( |14[ ) is of great importance in the remainder of 
this work, we now discuss the steps required for obtaining 
it. Starting from the quadratic Hamiltonian of Eq. (13) 
one first performs the Fourier transformation 


7V-1 


Vn ^ 


• 2 - 7 rjfc 




(15) 


k=0 


of the mode operators Ofc. After substituting Eq. (15) 


into Eq. (13) the Hamiltonian can be written in matrix 
form as 


g (.1, (4 J 




(16) 

This block-diagonal structure of the Hamiltonian in 
terms of the mode operators bk and bN-k makes it partic¬ 
ularly useful for calculating the operator O of Eq. (10). 


Finally, from Eq. (16) one can go to the diagonalized 
Hamiltonian of Eq. (141 by performing a suitable Bo- 


goliubov transformation on the two-dimensional block of 
mode operators bj., &Ar_fc 


Cfc = cos 9k bk - sin 9k 
CN-k = sin 9k bl + cos 9k b^-k 


(17) 


with Ok = ^ tan 


-1 




2 I Jcos(2^)+B I '^•5 2^ 


and (j)j = f ,Vj. 


We now determine the optimal achievable precision in 
estimating J given that B is known. 


with 


Oj{J,B,t)= [ U{J,B,s)HiU\j,B,s)ds, (18) 

Jo 


where Ffi — Writing the latter in terms 

of the fermionic operators bk one obtains 


H 


1 — 


N-l 

E 


[K bN-k) (^_,sin(^) -cos(^)J [b^j,_ 


Mu 


(19) 

In order to calculate the operator Oj{J,B,t) of 
Eq. (18) we need to determine the action of H{J,B) on 


the fermionic operators bk, i.e., we need to determine 
bk{s) = U{J,B,s)bk{0)U^J,B,s). This is simply the 
Heisenberg equation of motion for the mode operators 
bk- Due to the block-diagonal structure of H{J,B), when 
written in terms of fermionic operators bk, the solution 
to the Heisenberg equation of motion can be seen to be 


fefe(s) \ _ (_■ ( iPk\\ ( bk{0) 

&]v-fc(^)j" H '"[-^Pk -ak))\bl_ki0))- 

’V*' 

Rk{s) 

( 20 ) 

Henceforth, we drop the explicit time dependence of the 
mode operators bk for convenience. 


Substituting Eq. (201 into Eq. (18) yields 


Oj{J,B,t) = Y, (4 bN-k) ( f dsRl{s)MkRkis)) (J 

k=o ^ ^ ^ 


h 

'N-k 


Performing the integration over s gives 


( 21 ) 


A. Estimating the Interaction strength 

We now proceed to estimate the interaction strength 
J for the Ising Hamiltonian. The QFl is given by Eq. © 


'’"T (a) U\) 


where 


( 22 ) 


— 




Ok (ofe cos(^) + (3kSm{^))t 


gTTfc'i'i k _i_ 9k{Pk cos(^)-au sin(gf^)) mn{2ujut) 
2uJk 




ih {ak cos (^) + /3fc sin (^)) t + {h cos (^) - sin (^)) (sin^ (w^f) - 


(23) 


where we have separated the linear and oscillatory be- haviour of Vtk and (see also m) and uJk = y/cti+ Pi- 
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As Eq. (22) is of the same form as Eq. (16) it can be 


brought to the diagonal form 


N-l - 






(24) 


by a suitable Bogoliubov transformation (see Eq. (|1 
Note that the free fermions corresponding to 
different from those of Eq. (14), and we may choose with- 


d\, dk are 


out loss of generality the positive square square root of 
ill + which corresponds to choosing the vacuum 

state for free fermions to be zero. 

The optimal precision in estimating the interaction 
strength J is now easy to determine. As the latter is 
inversely proportional to the square root of the variance 


of Oj{J, B,t), we simply need to determine the maxi¬ 
mum achievable variance for the operator in Eq. (24). 


This is achieved by preparing the equally weighted su¬ 
perposition of the vacuum state and the state where all 
N modes are occupied by fermions. The variance with 
respect to this state is simply given by 




A'^Oj{J, B, t)max = 




|A, 


(25) 


. fc =0 


As the sum includes N summands, the variance of 
0,j{J, B, t) scales as up to some constant factor that 
depends solely on the ratio between J and B as we now 
explain. 


Using Eq. (23) one can easily show that 


^1 + 


AJ" = 


(/3fc cos (^) - ak sin (^))^ (1 - cos (2a;fet)) -h (ofc cos (^) -F jdk sin (^))" 


(26) 


For long interaction times, i.e., t —>■ oo the term quadratic 
in t in Eq. (26) completely dominates. Moreover, for N 

can be replaced, to a good 


(25) 


large the sum in Eq. 
approximation, by an integral resulting in the following 
simple expression for the variance 




G{g) = 


/ (f-HffCOs(a:))^ 
27r io V 1 + ^5 cos(a;) 


dx 


(27) 


Thus, the variance scales as i.e., at the Heisenberg 

limit, up to an overall constant factor that only depends 
on the ratio B/J. The function G is plotted in Fig. 
Notice that G{B / J) exhibits a phase transition at J = B. 

We now proceed to determine the optimal precision 
with which one can estimate J, knowing B, if we restrict 
the input state of N systems to be a product state. As it 
is not immediately evident what separable states look like 
in the basis that diagonalizes the operator Oj{J, B, t), we 
work directly with Eq. (fO). However, due to the form 


of Eq. (f 0) it is highly non-trivial to perform an analyt¬ 


ical optimization over all possible separable states of N 
qubits. To that end we perform a brute-force optimiza¬ 


tion of the variance of Eq. (25) over all possible input 


product states for up to A = 11 qubits. The results are 
shown in Fig. where one can already see a difference in 
scaling between the optimal product state strategy and 
the optimal quantum strategy even for small system sizes. 
Moreover, within a high margin of certainty the scaling 
of the QFI using the optimal product input state is at 
most linear in N. 



FIG. 1. The asymptotic value (large t and large N) of the 
normalized variance A^Oi/N'^B where i G (J, B), as the fnnc- 
tion of the parameter ratio g. The blue line corresponds to 
the variance of A‘^Oi/{Nt)^ with respect to the optimal state 
whereas the orange curve corresponds to the constant factor 
F{g) that multiplies the quadratic term of the variance of Oi 
with respect to the GHZ state. Observe the phase transition 
at the point where the parameters are equal. 


B. Estimating the field strength 

We now proceed to estimate the magnetic field B, 
given we kn ow J exactly. The procedure is identical to 


that of Sec. 


HI A 


Writing H 2 = 


the fermionic operators bk we obtain 


N-l 


■^2 - X! 


k=0 


1 0 
0 -1 


h 


’N-k 


in terms of 


(28) 
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Chissicci] vs Optiiiiril estimation of ./ 



FIG. 2. Comparison between A^Oj{J, (blue cir¬ 

cles) and the variance of Oj{J,B,t)/t^ using product input 
states (red squares). Both expressions are computed for the 
case J = B = 1 and t = 20. The blue circles represent the 
exact analytical value for A^Oj{J, B, t)max/t^ using Eq. (231. 
The minimum squared error fit for the red squares is given 
by Oj{J, B, = aN^ + c with a = 0.7476 ± 0.2815, b = 
1.034 ± 0.1350, c = -0.5139 ± 0.5831 with a 95% confidence. 
The green line represents the optimal QFI for the case where 
H = i.e., when B in Eq. Q is set to zero, 

using the optimal product state \ip) = |01)®^^. 


Substituting Eqs. (20 28) into Eq. (101 and performing 
the integration over s yields 

N-l 

^ (bl bpj-k 


k=0 


Ak Bk 
Bl -Ak 


bk 


b{ 


N-k 


where 


Ah = 


Bh — 


(29) 


alt Pi sin(2a;fet) 


UJt 




iakPkt , 2/3fca;fcSin^(wfet) - iafc/3fcsin(2a;fct) 






(30) 


where we have again separated the linear and oscillatory 
parts in t. As Eq. (29) is of the same form as Eq. (16) it 


can be brought to the diagonal form 

N-l 


OBiJ,B,t) = 2Y,\/Al + \Bk\^ flfk 


(31) 


A;=0 


by a suitable Bogoliubov transformation (see Eq. 0). 
The optimal achievable precision for estimating the field 
strength B is given by the optimal variance of Ob{J, B, t) 
in Eq. (31) which can be easily computed to be 

A^OB{J,B,tU,^= (f^^Al + \BkA , (32) 

\fc=o / 


and is again achieved by the equally weighted superpo¬ 
sition of the vacuum state and the state where all N 
fermionic modes are occupied. We note that because the 
Bogoliubov transformation diagonalizing Eq. (29) explic¬ 


itly depends on the coefficients of Eq. (30) the fermions 


described by modes fk here are different than those de¬ 
scribed by modes and in Eqs. (24 14) respectively. 
As the sum in Eq. 


includes N summands, the max¬ 
imum variance of Ob{J, B,t) scales as N'^ up to some 
factor. 

Using Eq. ([30| one can easily show that 


Ak 


\BhP = 


2 _ Pi (1 - COS {2uJkt)) + 2t^alujl 

2^^fc ■ i i 


For t —)■ 00 and N very large the quadratic term in t dom¬ 
inates and the maximal variance can be given explicitly 
as 


A^Ob{J, B, t)^ax = NH'^G (J/B) (34) 

where G is the function given in Eq. ( |27| ) and Fig. 

That the optimal precision in estimating either B or 
J is asymptotically given by the same expression can be 
understood via the duality of the one-dimensional Ising 
chain [33]. The duality is associated with whether one 
adopts the spin degrees of freedom on the chain or the 
kink degrees of freedom—associated with the links of the 
chain—as qubits. If adjacent spins are parallel then there 
is no kink, else there is a kink. One can then recast the 
Ising Hamiltonian of Eq. ([^ in terms of the kink degrees 
of freedom as 

N N-l 

H{J, H) = J ^ ri*) + B ^ , (35) 

i=l i-1 


where {tq : a € (x,y,z)} have the same commutation 
relations as the Pauli matrices for the spin degree of free¬ 
dom. Up to a basis change, the Hamiltonian in Eq. (35) is 
identical to that of Eq. ([^, except that the roles of J and 
B are reversed. This is the reason why for the estima¬ 
tion of J the constant factor is given by G{B/ J) whereas 
for the estimation of the field it is given by G{J/B). We 
note that the duality works well for spins and kinks in the 
bulk of the one-dimensional chain but is problematic with 
spins and kinks close to the edges of the chain. However, 
for large enough N the effects at the boundaries of the 
chain can be safely neglected. Also note that the quater- 
nionic representation of the operators in the kink degrees 
of freedom is different than that of spins in one crucial 
way. The degeneracies of eigenstates in the kink repre¬ 
sentation are not the same as the ones for spins. Indeed, 
one can easily show that the ordered spectrum of eigen¬ 


values of 

with eigenvalue m having a degeneracy of 
the ordered spectrum of i® given by {Am = 

^ — {m + ^)\m G {0,..., N — 1)} with corresponding 
degeneracies given by 2('^^^). Notice however that the 
duality is not exact for finite N. 


di) 


is (Am = - m| m e (0,..., A^)} 

whereas 


/N\ 
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This difference in the number of degenerate states can 
be exploited, in the case of finite TV, to improve the pre¬ 
cision with which one can estimate parameters in the 
Bayesian estimation regime. 

The optimal precision for estimating B, given we know 
J, using a separable strategy is again numerically calcu¬ 
lated for up to = 11 qubits with the results shown in 
Fig-i Just as in the case of estimating J, one can al- 


Cla^.sical vs Optimal estimation of D 



FIG. 3. Comparison between A^Ob{J, (blue cir¬ 
cles) and the optimal variance of Ob{J, using product 

input states (red squares). Both expressions are computed for 
the case J = B = 1 and t = 20. The blue circles represen t th e 
exact analytical value for A^Ob{J, B, t)max/t^ using Eq. (301. 
The minimum squared error fit for the red squares is given 
by Ob{J, B ,t)/B — aN^ -|- c with a = 1.099 ± 0.1350, b = 
0.9114 ±0.0423, c = -0.2747±0.2421 with a 95% confidence. 
The green line represents the optimal QFI for the case where 
H — B i.e., when J in Eq. Q is set to zero, using 

the optimal product state \ip) = |±)^ 


ready observe a difference in scaling of the QFI between 
the optimal product state strategy and the corresponding 
optimal quantum strategy. Moreover, within a high mar¬ 
gin of certainty, the scaling of the QFI using the optimal 
product state is at most linear in N. 


IV. PRECISION USING GHZ TYPE STATES 

In the previous section we showed that the optimal 
precision in estimating either the interaction strength 
or magnetic field is achieved by states that are linear 
superpositions of the vacuum and all N sites occupied 
by fermions. Whereas such states can be prepared ef¬ 
ficiently, i.e., with a quantum circuit that grown poly- 
nomially with the number of qubits |45j . preparing such 
states in practice may still be quite challenging due to the 
optimal states dependence on both time, and J (B). An 
important question, then, is whether there exist states 
that are easy to prepare and yield Heisenberg scaling in 


precision for all time and all values of J (B). For ex¬ 
ample in ion-trap set-ups one can prepare the GHZ state 
using a single Sprensen-Mplmer gate [15] ■ In this section 
we analytically determine the performance of GHZ-type 
states for estimating either the interaction strength or 
magnetic held. We will show that the GHZ-type states 
allow for Heisenberg scaling in precision for all values of 
J (B), with only a constant factor difference from the 
optimal precision achievable. 

Let us hrst determine the ultimate precision of esti¬ 
mating B using the GHZ state 

\GHZ) = ^ ± |1)®^) . (36) 

The GHZ state yields the ultimate precision in estimat¬ 
ing B for the case where the latter is imprinted via the 
unitary operator U = exp{—iBtH 2 ). In order to calcu¬ 
late its performance for the Ising Hamiltonian of Eq. ([^ 
we will use the Jordan- Wigner transformation and the 
discrete Fourier transform to express the GHZ state in 
terms of the fermionic creation operators b\.. 

We begin by noting that the GHZ state can be written 
in terms of the fermionic operators {a, al} as 




N 


(37) 


Using Eq. (151 to transform the fermionic operators, {a^} 
to and noting that the Fourier transform leaves the 
vacuum invariant, yields 


|Gi7Z) = ^(|0)±e*^|l)), (38) 

where 0 is a relative phase that arises due to the Fourier 
transform. 

Using Eq. (381 to calculate the variance of the operator 
in Eq. (291, one finds that 

{GHZ\ Ob{J, B, t) \GHZ) = 0 


fN-l 


{GHZ\ 0|( J, B, t) \GHZ) = ^ A, 


(39) 


, fc =0 


and consequently A‘^Ob{J, B,t) = 

{GHZ\Ol{J,B,t)\GHZ). Using Eq. (|^, denot¬ 
ing 6 = and letting V —)• oo the variance of 

A^Ob(J, i?,t) explicitly reads 


A^OB{J,B,t) = — 


(1 ± geos 9)'^ 


d0 


47’’^ \Jo I + cos 0 

' (40) 

where g = J/B. From Eq. (40) it follows that the vari¬ 
ance of Ob{J, B, t) with respect to the GHZ state always 
scales quadratically with N, i.e., at the Heisenberg limit, 
with a prefactor F{g) slightly lower than that for the 
optimal states as shown in Fig. 










We now determine the precision with which one can 
estimate J using the state 


100 )^ 


| 01 )' 


V2 


(41) 


which is, up to a basis change, the state that yields the 
ultimate precision in estimating J for the case where 
the latter is imprinted via the unitary dynamics U = 
exp(—zJ_ffi). We note that since Hi has a doubly de¬ 
generate eigenspace for both its minimum and maxi¬ 
mum eigenvalue, the optimal state for estimating J is 
not unique. 

If we express this state in terms of the kink degrees of 
freedom, as opposed to the spin degrees of freedom then 
what we obtain is, up to an overall Hadamard transfor¬ 
mation, the GHZ state, i.e, a linear superposition of zero 
kinks, and N kinks. As the Ising Hamiltonian expressed 
in the kink degrees of freedom is given by Eq. ( [3^ , up 
to an overall Hadamard transformation, it follows that 
the variance of Oj{J,B,t) with respect to the state \ip) 
is given by Eq. (40), with g = B/J. 


V. CONCLUSION 

In this work we investigated precision limits for noise¬ 
less quantum metrology in the presence of parameter de¬ 
pendent Hamiltonians, and in particular the Ising Hamil¬ 
tonian. We showed that the ultimate limit in estimating 
the interaction or magnetic field strength scales quadrat- 
ically with the number of probe systems N, i.e., at the 
Heisenberg limit, and that the states that achieve this 
precision are linear superposition of the vacuum and fully 


occupied states of N free fermions. Moreover, whereas 
the Ising Hamiltonian generates entanglement this entan¬ 
glement does not help in boosting the precision scaling 
with respect to N that can be achieved with product 
states. In addition, we showed that the achievable preci¬ 
sion in estimating either J or B for the Ising Hamiltonian 
using GHZ-type states also scales at the Heisenberg limit. 

Whereas we have shown that the entanglement gener¬ 
ating properties of the Ising Hamiltonian do not boost 
the precision in estimation for product states, it may 
still be the case that we can exploit this property of 
the Ising Hamiltonian to reduce the amount of entangle¬ 
ment required in the initial input state of the N probes. 
This would be of high interest for practical realizations 
of quantum metrology where the creation of highly en¬ 
tangled states of N systems remains a challenge. 

In addition, our analysis deals with optimal states and 
bounds in the absence of noise. It would be interest¬ 
ing to investigate the achievable precision bounds in the 
presence of several physical noise models, such as uncor¬ 
related dephasing or depolarizing noise, as well as spatial 
and temporal correlated noise. Furthermore, it would be 
interesting to determine which of these types of noise can 
we readily combat via the use of error-correcting tech¬ 
niques, or by dynamical decoupling [471-149]. 
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